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Abstract 

Resonant tunnelling is studied numerically and analytically with the 
help of a three-well quantum one-dimensional time-independent model. 
The simplest cases are considered where the three-well potential is poly- 
nomial or piecewise constant. 

1 Introduction 

Even though the expression "tunnel effect" was coined for the first time in 
1931 by Schottky in German, ( "Wellenmechanische Tunneleffekt" , according to 
Merzbacher, 2002), the importance of the quantum transmission through a po- 
tential barrier at an energy below its maximum was acknowledged immediately 
after quantum mechanics reached its maturity in 1926. Indeed, Hund published 
the first papers on tunnelling in 1927 where a one-dimensional double-well was 
introduced to modelise chemical binding (Hund, 1927a) and the deformation 
potential of a NHa-like molecule (Hund, 1927b). In the latter article, Hund ob- 
tained the crucial point that tunnelling was exponentially sensitive to the char- 
acteristic parameters of the barrier. The very same year, tunnelling was first 
considered in unbounded models by Nordheim (1927) to describe the electronic 
emission from metals. One year later, Fowler & Nordheim (1928) successfully 
described how this emission was driven by high electric field and Gamow (1928), 
Gurney & Condon (1928) showed how tunnelling was involved in the a-decay of 
some radioactive nuclei; in both cases, the exponential dependence of the decay 
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rates with energy was explained by tunnelling through a one-dimensional bar- 
rier that modelled the mean attractive field created respectively by the metal 
and the nucleons (and the latters were still unknown in 1928). 

Such tunnelling models can be introduced very early in a first course on quan- 
tum mechanics. As soon as the stationnary Schrodinger equation is presented 
together with the constraints imposed on its resolution (normalisation, smooth- 
ness of the solutions), tunnelling can be understood as a direct consequence of 
the wavy nature of quantum particles. For usual textbooks treatments see for in- 
stance (Bohm, 1951, § 11.5), (Messiah, 1959, chap. Ill, § 7), (Merzbacher, 1970, 
chaps. 5 and 6). Probably more striking for the beginners than the quantization 
of the energy levels of a bounded system, tunnelling provides a first contact 
with the strangeness of quantum phenomena, even for those already familiar 
with evanescent waves in optics. 

Tunnelling has nowadays been extended to qualify any quantum process 
that is forbidden by real solutions of classical equation.^ thus, since a quan- 
tum/classical comparison is somehow necessary, it is, by definition, a phe- 
nomenon that occurs in the semiclassical limit, where formally fi. — > (physically 
this corresponds to large classical actions compared to the Planck constant. This 
limit is not trivial : from the Schrodinger equation we can see that the value 
fi. = is singular; however, in practice, the semiclassical approximation works 
quite well even for relatively large value of H. The only requirement in our case 
is that h must remain small enough to get two almost degenerate eigenvalues 
below the tunnelling barrier.) With the help of approximate methods of resolu- 
tion of ordinary differential equation that were originally developed by Jeffreys 
(1925), Kramers (1926), Brillouin (1926), Wentzel (1926)— the so-called JWKB 
approximation that provide h asymptotic expansions for the quantum quantities 
(for a rigourous though accessible presentation with more details than the ones 
usually presented in quantum mechanics textbooks see Bender & Orszag, 1978, 
chap. 10) — one may also gain some insight of how subtle the transition be- 
tween the classical and the quantum world occurs. However, from the singular 
behaviour of the wave functions (in their semiclassical amplitudes near the clas- 
sical turning points and, in any case, in their phases), it appears immediately 
that the semiclassical limit ?i — >■ is not as simple as the other familiar limit 
c — > CXI that describes the transition towards non-relativistic theories (physically 
this corresponds to small classical velocities compared to the speed of light c) . 

Remaining at an introductory level, within JWKB approximation, it is quite 
simple to justify the exponential form taken by the transmission factor t (de- 
fined to be the square of the ratio between the amplitudes of the outgoing and 
incoming waves) of one barrier like the one depicted in figure [l}i) (Bohm, 1951, 
§ 12.13), (Landau & Lifshitz, 1958, § 50, eq. (50.5)), (Messiah, 1959, chap. VI, 
§ 10 and exercise 5) (Merzbacher, 1970, chap. 7, eq. (7.33)) 

r ^ ae-2^/\ (1) 

or the tunnelling beating period T = h/AE in a symmetric double-well like the 
one given in figure [Dd), (Landau & Lifshitz, 1958, § 50, problem 3). In that 

^For more than one degree of freedom or for time-dependent systems, the interdiction may 
come from constraints other than the energy conservation. If we look at the dynamics in the 
phase-space and make the canonical transformation that exchanges position with momenta, 
it also appears that the reflexion above a barrier can be considered as a tunnelling process. 
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Figure 1: a) One simple potential barrier and b) the double- well potential pro- 
vide the first illustrations of tunnelling in one dimension. 



latter case, tunnelling appears as a coupling between two states localised in the 
two wells that lifts the degeneracy of their energy. Actually, below the maximum 
of the barrier T^ax the symmetric and antisymmetric exact eigenvalues come 
by pairs and E~ labelled by a finite number of integers n = 0,1,..., N and 
the splitting of the nth doublet is given by 

AEr^^^'lE^-E-] ahe-^/\ (2) 

It requires some skill to get the correct prefactors a (Garg, 2000) but if we do not 
bother about them, the above expressions can be straightforwardly obtained. 
In both cases A is the classical action "under the barrier V{q)" of the particle 
with mass m, more precisely the imaginary part of the action computed with 
one branch of the imaginary momentum: 

A{E) II' / j2m{V{q) - E) dq , (3) 

where qr and q'^ are the positions of the turning points at the energy E at 
which the barrier penetration occurs. In the double well case, E is the average 
energy of the tunnelling doublet En = {E^ + E~)/2 and may be determined 
semiclassically within one well by the Einstein-Brillouin-Keller quantization con- 
dition (Bohm, 1951, § 11.5), (Landau & Lifshitz, 1958, § 48), (Messiah, 1959, 
chap. VI,§ 11), (Merzbacher, 1970, chap. 7,§ 3), 

2nh(n+\) (4) 

r<lr(E) . 

2/ J2m{E-V{q))dq, (5) 

where q'^{E) stands for the most right turning point at energy E. 



p{q) dq = 

closed trajectory 
at energy E 
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Simple exponential laws like ((J) or discovered by the pioneers of tun- 
nelling mentioned above, are advantageously used in many applications, most 
notably the principle of the scanning tunnelling microscope (or rather nanoscope) 
(Eigler & Schweizer, 1990). However, when the potential has a richer structure 
than just one bump, even when staying in one dimension, strong deviations 
by several order of magnitude on the transmission or on the splitting are to be 
expected. For instance, even if the transmission of one barrier can be made arbi- 
trarily small because of ([T]), the total transmission factor of two such successive 
bumps may be enhanced up to its maximum value 1. This resonance tunnelling 
has fundamentally the same roots than the resonance scattering by one well 
(Bohm, 1951, chap. 11, § 7), (Messiah, 1959, chap. Ill, § 6) (Merzbacher, 1970, 
chap. VI, § 8): an enhancement is produced by the constructive interference 
of waves reflected back and forth in between the two bumps. In optics such a 
coherent superposition of "trapped" electromagnetic waves is the principle of 
the Fabry-Perot interferometer (Born & Wolf, 1980, § 7.6). Though accessible 
with no more technical nor conceptual tools than the ones involved in the jwkb 
derivation of (P), ([2]) and quantum analogous treatment of the Fabry-Perot 
effect is scarcely treated in the introductory literature (Bohm, 1951, chap. 12, 
§§ 14-17), see Granot's (2006) recent publication in this journal. The aim of 
the present article is to fill a gap and complete the study of resonant tunnelling 
by considering the spectrum of bounded systems rather than scattering. Mul- 
tidimensionnal tunnelling (including time-dependent one-dimensional systems) 
requires theoretical tools far beyond the readership's interest of the present jour- 
nal and is still the subject of an intense research field. However, recent studies 
has shown that resonant tunnelling (mainly in bounded systems) paved the 
road to a better understanding of tunnelling in complex systems (Schlagheck 
et al, 2011). For instance, the huge fiuctuations by several orders of magnitude 
in some tunnelling doublets are precisely of the same origin than the resonance 
that will be carefully described in the present article. We will remain with a 
one-dimensional static potential V and choose two tractable models that in- 
volve three wells. First (section [3]) we will consider a smooth potential, the 
simplest one being given by a polynomial of degree 6. One of its advantage 
is that the spectrum can be numerically computed with very fine precision at 
low cost. Then (section [3]) we will introduce another model where y is a piece- 
wise constant which will allow us to go further in the analytical computations. 
Comparison between the two will strengthen the interpretation of resonnant 
tunnelling in terms of avoided degeneracies in the spectrum and will show that 
the Fabry-Perot effect provides a correct interpretation. The two models en- 
capsulate the same physics and, indeed, their comparison may help to support 
this intuition. What seems more interesting for pedagogical sake is to con- 
sider the opposite way : starting with the intuition that the two models are 
physically qualitatively equivalent, the very different methods to study them 
(numerics, level dynamics, Husimi distribution for the smooth case and analyt- 
ical computations, semiclassics, transfer matrices for the other one) appear to 
be complementary and mutually supporting. 

All along this paper, we will use mainly physical concepts and mathematical 
tools that correspond to a undergraduate course on quantum physics. However, 
the accumulation of different techniques (specially the somehow lengthy compu- 
tations of section 3.2, a reasonable skill in asymptotic expansions, the numerical 
implementation and the notion of quasiclassical states) are more accessible once 
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the very first notions of quantum physics are well-understood or during a special 
supervised training course. 



2 Three-well polynomial potential 
2.1 Presentation 

The simplest smooth symmetric one-dimensional bounded potential with reso- 
nance is given by 

V{q) = {q^^a'f{q^-h'), (6) 

where a and h are real. For a suitable choice of the parameters (a, &), with 
& > a, the potential exhibits two symmetric external wells separated by a third 
deeper one (Fig. [5^). In the following, we will systematically work with a non 
relativistic particle whose mass is taken to unity. The quantum spectrum is 
purely discrete, bounded from below and the eigenstate \'ip) of energy E is given 
by the spectral equation 

^IV') = E\i,) (7) 

where the Hamiltonian is if = H{p^ q) where p, q are the momentum and posi- 
tion operators and H the classical Hamiltonian 

H{p,q)^^ + V{q). (8) 



2.2 The spectrum 



When dealing with a polynomial Hamiltonian, a convenient way to perform 
the numerical computations is to write directly the canonical operators in the 
harmonic oscillator eigenbasis (Korsch & Gliick, 2002) 
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_ (9) 

The Hamiltonian ([8]) is now written in terms of power of these infinite matrices. 
In practice, one needs to truncate them. First, one computes numerically the 
products of the truncated matrices and then one sums each terms. One can 
now diagonalize the Hamiltonian to get the eigenvalues. Numerically, we must 
finally check that the truncation of these infinite matrices does not affect the 
results up to the required precision. 

The eigenenergies are plotted as a function of the Planck's constant h in 
figure ^Bp) . Even though it is a standard procedure when one studies quan- 
tum phenomena in the semiclassical limit, it is worth to remind that treating 
the Planck constant as a varying parameter has to be understood as a substi- 
tute for varying the physical parameters of the system (for instance, through 
a rescaling of the quantities including h in order to work with dimensionless 
parameters). Since V{—q) = V{q), the eigenstates can be classified according 
to their parity (±). Moreover, for small enough h, the part of the spectrum 
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below the local maximum Vmax can be divided in two families: first, a family 
of doublets , E~ ) that remain positive and below Knax for any value of h 
(like the red lines in fig. [Sb)) and second, a family of energies that eventually 
become negative for small enough H (like the blue line in fig. [Dd)). Except near 
the special values of H where the two families almost intersect (these crossings 
will be discussed in the next section) the first family corresponds to states lo- 
calised in the lateral wells with ~ ^ huji{n + 1/2), with integer n 
and uji = y^V"{a) the harmonic frequency in the bottom of the lateral wells : 

like in the double- well case, the difference Ai?„ =^ — E^\ will be our tun- 
nelling splittings between symmetric and antisymmetric states. The energies 
that make up the second family do not form doublets but rather correspond 
to states localised in the central well. In the semiclassical limit these energies 
behave as E^^^ = V{0) + hLJc{m + l/2) + O{h'^) with integer m and lOc = ^^"(0) 
the harmonic frequency in the bottom of the central well. 




Figure 2: a) The three-well potential ^ for a = 1.5 and b = 0.6. b) The first 
20 energy levels versus 1/h (black lines in figure b). The dashed green lines 
delimit the energy range where the doublets Ai?„ lie. We have also sketched 
the doublet AEq (in red) and the level iJg^'' (in blue) associated to the central 
well. The inset in b) shows the avoided crossing between the three levels E^ , 
E^ and £'5'^'' that occurs near l/?i ~ 8.18. 



2.3 Avoided crossings in the spectra 

In the absence of singularities (Cohen & Kuharetz, 1993), the spectrum has no 
degeneracies. When varying a control parameter {h in our case), what appears 
to be a crossing between energies is actually an avoided crossing if one look 
at a sufficiently high resolution (see the insert in figure [5Jd)). More precisely, 
two energy levels having the same parity seem to repel one from the other and 
a substitution between one component of the doublet (Eq) and E^'^^ occurs 
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when 1/fi increases from 8.14 to 8.24. This scenario can be reproduced using a 
three -level model f? l3.3p and can be followed using a phase-space representation 
of the corresponding wave functions (ij 12.41) . 

Let us note that near such avoided crossings, the definition of the doublet 
becomes ambiguous: whatever criterion we retain to discriminate between the 
levels and Em^ (for instance the two nearest energies or those states selected 
from their overlap with a coherent state localised in one of the three wells) a 
discontinuous switch will occur and the level repulsion originates a large in- 
crease, by several of magnitude, of the tunnelling splittings Ai?„. Those are the 
resonances seen from the energy levels point of view. The word "resonance" is 
then to be taken in the usual sense: the strong variation of AE within a small 
range of variation of one parameter (l/h) is directly related to the coincidence 
between frequencies, namely the frequencies of the Rabi oscillations between 
the wells. From then on, we will define AEn in selecting the energies involved 
as the largest overlap between the eigenstates and a coherent state localised in 
one of the external well. The splitting will thus be computed numerically fol- 
lowing this criterion, leading to the figures |3^) andlHl As seen in figure |3^), the 
simple exponential behaviour ^ breaks down when the doublet is mixed with 
a third level and some resonant spikes appear. Moreover, these peaks are not 
differentiable at the values of h corresponding to the discontinuous switches. As 
AE increases close to a resonance, the tunnelling period T — h/AE becomes 
smaller: a quantum state initially prepared in one of the external wells will 
then oscillate more rapidly from this well to its symmetric. For some particular 
values of the parameters, that is to say when an isolated level comes to perturb 
the doublet, a middle well separating two symmetric wells can thus enhanced 
significantly the tunnelling. In section [3] we will derivate almost analytical 
formulas for the splitting and the height of the peaks. 



2.4 Husimi representation 

It is very illuminating to see the resonance at work using a phase-space repre- 
sentation of the eigenfunctions involved in the triplet E^ , E^n' near an avoided 
crossing. Instead of computing eigenstates in either the q— or p— representation, 
we will choose the Husimi distribution ipiP^ q) which mixes both space and mo- 
mentum variables and gives the density probability to find a quantum state in 
the phase space. The Husimi distribution is defined as the square modulus of 
the projection of the state onto a coherent state \a) (Cohen- Tannoudji, Diu 
& Laloe, 1997, complement Gy). One can use the basis \n) of the harmonic 
oscillator in order to write the coherent state as follows 



- P-|"IV2 



E 

n=0 



(10) 



dcf 



where a — {p + iq)/\/2h, then 



'>Pip,q) =^ \ {a\ip)\ 



2 



n=0 



(11) 



•^In more complicated systems which are beyond the scope of this paper, one can observe 
crossings in the spectra and then the sphtting would exhibit antipeaks with infinite height in 
a semilogarithmic scale. The tunnelling would thus be completely destroyed. 
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Figure 3: The splitting AEq is plotted in a) in a semilogarithmic scale versus 1 /H 
for the three-well potential (jH]) with a — 1.5 and b = 0.6. The figure b) exhibits 
the associated spectra versus 1/h. The double arrowed red line indicates the 
correspondance between avoided crossing in the spectra and resonances in the 
splitting. The resonance at ?i ~ 1/8.18 has been magnified in the insert in figure 
a). 



The Husimi distribution ip{p,q) is definite positive. The coherent state \a) is 
a Gaussian wavepackct built such that it saturates the uncertainty inequalities 
i.e. ApAq = h/2. A coherent state, also called a quasi-classical state, can be 
understood as the most classical quantum state since it can be shown that the 
expectation values of the quantum observables follow the same evolution as the 
corresponding classical ones at the leading order of h. Physically, the square of 
the scalar product of the eigenstate \4>) and a coherent state can be interpreted 
as the density probability to find the quantum state in a cell of area H centered 
at the momentum p and the position g in a coarse-grained phase space. The 
Husimi distribution appears thus to be an appropriate picture to study quantum 
dynamics in the semiclassical regime (Novaes, 2003). Far from a resonance (see 
figureS^.O), the symmetric and antisymmetric eigenstates related to the doublet 
with the energies {Eq , Eq ) are mainly localised in the external wells while the 
isolated level is clearly well-delimited in the middle well as shown in figures 
11^.1,2,3). When 1/h increases from 8.15 to 8.25, we can follow the exchange 
between the state corresponding to -Eg*^"* and the state corresponding to E^^ while 
the even state remains almost unaltered and mainly localised in the lateral wells. 
At the resonance (fig. |3}d.1,2,3), the two states with the same parity become 
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Figure 4: For the three-well polynomial potential (IHl) with a = 1.5 and b = 0.6, 
we have plotted the Husimi distribution for the three levels involved in the 
avoided crossing around h ~ 1/8.18 (see also insert in figure [2^). The vertical 
blue lines in respectively a,b,c.O) defines the values of the Planck's constant 
(resp. l/h~ 8.14, 8.184, 8.24) and shows the relative position of the energies in 
the spectra (black lines) where the dashed arrows point out the energy levels 
related to each Husimi distribution depicted in a,b,c.l,2,3). The classical phase 
space is also superimposed on the Husimi plots. 

intertwined: this is the signature in the phase space of the so-called resonances. 
For these particular values of the parameters, the odd quantum states leak out 
through the barriers all over the three wells. We cannot clearly distinguish which 
odd state is now a part of the doublet and we find the ambiguity concerning the 



9 



definition of tlie splitting again. 



3 Almost analytical computations 
3.1 Square potential 

To go beyond a simple diagonalization and understand analytically the origin 
of the resonant spikes in the splittings, potential ^ is still difficult to handle. 
The standard instanton procedure of Wick rotating the time 1 1— )■ ii in order to 
work with a upside-down potential (Coleman, 1985, chap. 7) fails because there 
remain a barrier to cross when reversing V into — in figure[2]a). Two of us have 
recently shown how to generalise the instanton method in the most general Id- 
case (Le Deunff & Mouchet, 2010) and how to recover analytically the effect of 
resonances but this method goes far beyond the present elementary complement 
to a first years lecture on quantum physics. To deal with a simpler analytically 
tractable model let us introduce a three-well piecewise constant potential that 
will mimic the situation described in the previous section. Therefore we will 
take (Fig. [S]) 



-l-oo Region 1: q < —c, 

Region 2: — c < q < —6, 

Kiax Region 3: -b <q< -a, 

V{q) = { Klin Region 4: - a < q < a, (12) 

Vmax Region 5: a < q < b, 

Region 6: b < q < c, 

^+oo Region 7: c < q. 

While keeping the essential features of resonances described above, we will then 
be able to derive analytical formulas in the semiclassical limit for the splitting 
and the height of the resonance peaks. 
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Figure 5: The piecewise three- well square potential defined in (|12l) with the 
parameters a = 0.55, b = 1.25, c = 1.75, Kiax — 1-0 and Vmin ~ —1.82. 
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3.2 Semiclassical transfer matrix approach 

In the case of one-dimensional systems, a convenient way to obtain the eigen- 
functions and eigenenergies is to use transfer matrix formahsm (Walker & Gath- 
wright, 1994) (Azbel, 1983, Appendix A). The transfer matrices connect the 
amplitudes of the wavefunctions in the different regions of the potential. As 
we want to compute tunnelling splittings, we will pay particular attention to 
energies such that < E < V,nax- We will use this approach to get the equa- 
tions which give the energy levels and secondly to compute the splitting in the 
semiclassical regime. 



3.2.1 Energy quantization 



With potential (fT^ . the general solution ^p{q) of the time-independent Schrodinger's 



equation at energy E, —h'^ip"{q)/2 4 


- ViqMq) - 


: E'>jj{q), (the prime stands for 


the derivative) has the form 










'Region 1 




= 0, 






Region 2 




= yl2c"="9 + 


B2e-"=o9, 




Region 3 




= yl3c'^»"9 


+ ^36-'=— «, 


< 


Region 4 




= yl4e"=-"9 


+ S4e-''="'-«, (13) 




Region 5 


^b{q) 


= yl5e'=»-9 


+ B5e-''™9, 




Region 6 


ipeiq) 


= Aee''""' + 


B6e-"=°«, 




, Region 7 




= 0, 





where 



fco 



/2E 



y/2{E - Knin) 



\/2(Vniax — E) 



(14) 



The amplitudes (Ai, Bi) are a priori complex numbers. As the energy _E > 0, the 
wavefunctions ^"2(9), tpiiq) and ipe{q) are linear combinations of two plane waves 
travelling in opposite directions, while in the regions where the energy is lower 
than the barrier, ipsil) a^id ''Psil) are combinations of real exponentials. Though 
the potential V{q) is not continuous, one still requires the eigenfunctions and 
its derivative to be smooth everywhere and especially at each frontier between 
two regions, i.e. q = {±a, ±6, ±c}. It leads to the following equalities 



V'3(-«) = V'4(- 

i^',{b)=i^'e{b), 



b), 



(15) 



Using these continuity equations, one can express the amplitudes {Ai,Bi) in 
terms of any others {Aj , Bj ) through the relation 



'A, 

.B, 



T(*,j) 



T(*,j) 



Tii(i,j) 
T2l(^,J) 



Ti2(i,i) 
T22(*,i) 



(16) 
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where T(i, j) is the transfer matrix from the region j to i. Each transfer matrix 
can be written using three elementary matrices 




ika Q \ 

e--j' ""('^ 
1 /l-ifc/fc l + ifc/fc^ 
2ll + ifc/fc 1-ifc/fcj 



-ka 



(17) 



where fc, k and a are reaL The matrices u (resp. v) correspond to the propa- 
gation in a classicaly allowed region {E > V) (resp. forbidden region, E < V); 
the matrices w correspond to the continuity equations at one step. We have 

1 -1 



U„(fc)u^(fc) = Ua+p{k). 



(18) 



and w fulfills the condition [w{k,k)] ^ = w*{k,k), for real k and k, where w* 
denotes the complex conjugate of w. The transfer matrices are given by 



T(2,3) = Uf,(fco)w(fc 

max 1 max J ; 

T(3,4) = V_Q(fcmax)w*(fci„in,fcmax)u_a(fcn 
T(4, 5) = U_a(fcmin)w(fci„ax,fcmin)v-a(fcma 
T(5,6) = Vb(fcmax)w*(fco,fcmax)u6(fco)- 



(19) 
(20) 
(21) 
(22) 



Then, the relation between (A2,i?2) and (yl4,i?4) are given by the transfer 
matrix : 



T(2,4) = T(2,3)T(3,4) 

= U;,(fco)w(fc,„ax, fco)vb(fcmax) 



Ua ( fcmin ) W ( fcmax j ^min ) Va ( fcmax ) 



(23) 



A straightforward computation shows that the elements of the matrix trans- 
fer through a barrier have to be such that T22 = (Tn)* and T21 = (T12)*. 
Since the potential is symmetric, the matrix T(2,4) contains already all the 
information about the system. Indeed, noting that 



,4) = [T(4,6)]"1 = [T(5,6)]-MT(4,5)] 



(24) 



U_b(fco)v^f(fcmax, fco)V-f,(fcmax) 



U — a ( kmin ) ( fcmax j ^min ) V— a ( ^max ) 



and performing the symmetry-axis transformation (a, b, c) — > (—a, —b, — c), the 
expression is nothing but T(2,4). We distinguish the parity of the eigen- 
states with the condition ^^4(0) = for the even states while the odd ones fulfill 
?/'4(0) — 0. Using the continuity condition at q = — c, it leads to the algebraic 
system 



Tn(2,4) Ti2(2,4) 
(Ti2(2,4))* (Tn(2,4))' 



±^4 



(25) 



where ± stands for the parity of the eigenstate. The transcendantal equations 
for the energy levels are thus obtained easily by manipulating the system and 
eliminating the amplitudes from the previous equalities: 

Re [(Tii(2,4) -t-Ti2(2,4))e-"=«'=] = : even states, (26) 
Im [(Tii(2,4) - Ti2(2,4))e-"=«^] = : odd states. (27) 
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These conditions can be written in a more tractable form for the even states 

D+{E) = F+iE) + G+(£;)e-''=--(''-"' = 0, (28) 
and the odd states 

D_ (E) =^ F_ (E) - G- (£:)e-2'=m-(b-") = o, (29) 

with 

F+{E) =^ (ko cos[fco(c - b)] + fc„,ax sin[fco(c - b)]] (30) 



X ^fcmax COs(fcmin^j) ^min sin(fcffiinfl) 

G+{E) (^fcocos[fco(c- 6)] - fc„,axSin[fco(c- 6)]) (31) 

X ^fcmax COs(fcmin^j) + ^min sin(fcffii,iCl^ 

F_(£;) f fco cos[fco(c - b)] + fc,„ax sin[fco(c - b)]] (32) 



X ^fcmin C0s(fcniin'3^) + ^max sin(fcijiinfl) 

G-{E) f fco cos[fco(c - b)] - fc„,ax sin[fco(c - b)]) (33) 



X ^fcmin COs(fcniinfl) ^max sin(fcijiinfl) J • 

The zeroes of D± [E) give the exact discrete energies E — E^ for the three- weh 
square potential. 

In the semiclassical limit, the height of the barrier is much larger than the 
energies and the term of order one in the sums (^51) and ([^ dominate the 
decreasing exponential, so the energies can be approximated keeping only the 
first term 

D+{E) ^ F+{E)^Q, (34) 
D^E) ^ F^{E)^Q. (35) 

Within this approximation, tunnelling is neglected. There are two possible ways 
to cancel (l34|) or ((35)) as the functions F+ [E) and F- {E) are written as a product 
of two functions 

F+{E) = F,{E)xFr^,+ {E), (36) 
F^{E) = F,{E) X Fm^^iE), (37) 



with 



F,{E) =^ fco cos[fco(c - b)] + fcniax sin[fco(c - b)], 

Fm^+{E) =^ fcmaxCOs(fcmina) - fcminSin(fcmina), (38) 
Fm,-{E) = fcmin COs(fci„ina) + a). 
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The factor (E) appears in both cases (even and odd) and corresponds to 
the quantization condition in the lateral wells in the limit & — a — )■ +cxd and pro- 
vide the average position e„ of the doublets. The second factors Fm,±{E) in ([55]) 
and ([571) are nothing but the transcendantal equations which lead respectively 
to the even and odd energy levels of the finite square well (Merzbacher, 1970, 
chap. 6, § 8) and approximate the energies related to the central well of (IT^ . 
It is clear that, in the limit where we do not care about tunnelling, the energies 
associated to the external wells are degenerate as both quantization conditions 
and share the same first factor. Avoided crossings also disappear as the 
two factors of Fj^{E) (resp. F^{E)) in (p6)) (resp. (jSTI) ) vanish simultaneously 
for some particular valuesO of h. 



3.2.2 Semiclassical formulas for the splitting 

To recover the exponentially fine structure of the spectra that caracterise tun- 
nelling, let us expand the right hand side of the exact quantization conditions 
((28| and ((29|) around the average value of the doublet that has been determined 
above E^ — Cn + 

i?+(e„ + e+) ~ i?+(e„) + e+i?V(e„)+o((e+)2) =0, (39) 
I?_(e„ + e-) ~ i?_(e„) + e,T^-(e„)+o((e-)2) = 0, (40) 



Keeping only the exponentially dominant terms, we get 



Ai?„ = |i?--i?+| = |e--6+| ~ 



G_(e„) , G+(e„) 



Fi(e„) F|(e„) 



where Fj_(e„) — F^{en) x Fm,±{e„). One recovers on average the exponentially 
small behaviour of the splitting. Noting that p ~ hk where k is equal to the 
expressions in (1141) , the argument of the exponential is nothing but the classical 
action under the two barriers of the potential 



/—a t^h />b 

pdq+ pdq = 2 pdq = 2nk„,^^{b-a). (42) 
-b J a J a 

While F^{en) ^ 0, one denominator in (|^T|) may still vanish when Fm,±{en) = 
which provides a condition for determining the resonances. To fix the diver- 
gences in (j4ip and get a finite value of the height of the spikes, we have to go 
to the next order in the Taylor expansions (p9|) and (|40|. We distinguish both 
cases depending on whether the third level is even or odd. If the resonance is 
due to an even level e+ then i^^(e„) — and we need to expand I?+(e„ -I- e+) 
to the second order. The splitting is of order o[e~^'^'"'"'^''~''''] but it is no longer 
true close to a resonance where it becomes much larger because of the prefactor 
which is no longer of order one. The dominant terms are thus now of order 
o[(e+)^] and o[e~^'^"="'^''~''^] while it is still sufficient to keep only the first order 
term for the odd condition Z?_(e„ + e~), as previously. In case an odd energy 
level is responsible for a divergence in (j41[) . we proceed similarly with the odd 



■^For these particular values of h, some zeros of the conditions II34II and II35II are of order 2 
and give the average positions of the actual avoided crossings, otherwise they are all of order 



one. 
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condition D^{en + e„ ). Collecting all the results we get 



and 



G_(e„) 



-2fc„ax(f>-a) 



Fi(e„) 



^ -2G+(e„) 



-femax(b-a) 



^2G+(e„) 



(43) 



.(b-a) 



/2G_(e„) 



F:'(e„) 



K(b-a) 



(44) 



where AEn.+ (resp. Ai5^„^_) is the height of a resonance caused by an even 
(resp. odd) level from the central well. The denominators F!^{en) can be again 
simplified and written as the product 2Fj(e„) x _|_(e„) when the level is 
equal to e„ (at the resonance) . The argument of the exponential in the splitting 
formula (pi far from resonances is twice the argument in (|43l) and ([33]). The 
semiclassical formula (|^T|) shows a very good agreement with the exact results 
as plotted in figure |6l The semiclassical predictions are more and more accurate 
for large values of even around the resonances where the height of the peaks 
is very well predicted by the formulas (P5)) and (|44p . 



3.3 Matrix model 

We want to obtain a matrix model of tunnelling in a potential with several 
wells, as the one of equation ([B]). A solution is to represent the Hamiltonian in 
a basis of " quasi- modes" coupled by tunnelling. Clearly, tunnelling is a wave 
phenomenon which couples two adjacent wells classically separated by a finite 
barrier. If there was no tunnelling at all, then each well would have its own 
independant eigenmodes, approximately given by jwkb ansatz based on the 
real quantizing classical tori. In the semiclassical limit, tunnelling is a vanishing 
effect, and these jwkb modes become very good local approximations of the 
exact solutions of the time independant Schrodinger equation, that is, the exact 
eigenmodes. For this reason, these jwkb ansatz are often called "quasi-mode" in 
this context. But, if one sticks to real classical tori, then tunnelling is not taken 
into account by these approximations. On the other hand, one can consider 
tunnelling as a small coupling between quasi-modes. 

For each well of the potential, we hence define a quasi-modes family 0„ , with 
quasi-energies E^. In the simple case of a symmetric double- well, one has quasi- 
modes with energies i?^, associated with the left well; and quasi-modes (p"^ 
with energies E^, associated with the right well. Because of classical symmetry, 
both are degenerate, that is, i?^ = but tunnelling lifts this degeneracy in the 
exact spectrum. Close to energy E'^, the Hamiltonian is well described by a two 
level system in the basis of the two corresponding quasi-modes. For instance, 
close to the ground state, one has 
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Figure 6: The splitting for the ground state is plotted in a semilogarithmic 
scale versus \/h for the piecewise constant potential defined in with the 
parameters a = 0.55, b = 1.25, c = 1.75, Vmax — 1-0 and Vmin — —1.82. 
The black line is the exact computation obtained by solving numerically the 
exact quantization conditions ((28|) and (|29)) . The red line corresponds to the 
semiclassical formula (|4T]) and the green dots show the height of the resonances 
computed with (1^5]) and (|33]). The insert is a magnification of the resonance at 
n~ 1/9.093. 



where 5 gives the coupling strength induced by tunnelling between the two 
quasi-modes. Diagonalisation of this matrix gives two eigenstates, an odd one 
and an even one, 

V^± = -i= ± ^'D , (46) 
associated with two eigenenergies 

E±=Eq± 5. (47) 
In this model, the tunnel splitting is therefore given by 

A£: = 25. (48) 
Consistency with jwkb analysis and equation ^ is obtained by chosing 25 — 

It is easy to generalise this simple model to a triple well like ([6]) by defining 
three families of quasi-modes : (j)''^, (p^^ and cj)!^. Now we focus on the situation 
shown on figure SI In the language of this model, it is an avoided crossing 
involving a doublet (0jj,0^) with degenerate quasi-energy E^^, and a state 
with quasi-energy E!f^ which is close to El^. As H is tuned, E^ crosses El^, 
but tunnel coupling makes it an avoided-crossing. In the region of this avoided 
crossing, one can restrict the Hilbert space, in the same spirit as in (Tomsovic 
& UUmo, 1994), to a subspace spanned by these three states, that we now name 
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without label in order to simplify notation : (jji and (j)r with energy E, and (j)c 
with energy E' . In this basis one has something like 




~ (5 E' 6 \ . (49) 



Notice that there is no direct coupling from the left well to the right one, and 

only remains tunnelling from left to centre and from centre to right. In order 
to match with § 12.41 realistic values of these coefhcients would be 

E = hoji{n+^) (50) 

E' = V{0) + hoJc{m+^) (51) 

6 = ahe-"^/^. (52) 
It is easy to compute the corresponding eigenenergies : 

E^ = ^±^-l^iE~E'y+8S^ (53) 

E2 ^ E (54) 
E + E' 



E3 - —^ + ^ViE-Er+8d^. (55) 
In figure[71 we have plotted these values as a function of 1/fi, by using expressions 

dini), (EH) and dm). 

When E' - E^ 5, one has 

El ~ E - (56) 

2\E-E'\ ^ ' 

E2 = E (57) 

This corresponds to the situation of figure |3^.0) where only direct tunnelling is 
involved, and the central state has little influence. The splitting of the doublet 
E2 — El is of the order 5^, that is, e^"^^/^, which is basically the product of two 
transfer coefficients through one barrier. 
When E - E' ^ 6, one has 

El ~ £^'--r-- r (59) 

2\E-E'\ ^ ' 

E2 ^ E (60) 

This time, the energy of the central state (pc has crossed the doublet, which now 
corresponds to E^ — E2, as in the situation of figure S^. 3). 
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E 



0.305 



0.315 



0.31 




9.4 



9.42 



9.44 



9.46 



1/h 



Figure 7: Avoided crossing modelised with a 3 levels system. We have used 
uji = 5.833, ujc = 3.656, n = 0, m = 5, 1/(0) = -1.8225, a = 0.197 and A = 0.34. 
The vertical black dashed line indicates the value of h at the resonance condition 
E — E' . The numerical values are chosen according to the parameters used in 
figures [5] and HI 

The situation of figure |1^.2) corresponds to the value of h for which the 
central state has resonance E = E' with the lateral ones, and one then has 



The splitting then reaches its maximum l^E = 2^/25, of the order e"^/'* which 
is consistent with (|43p and (jH]). Resonance with the central state thus facilitates 
tunnelling, by one order of magnitude, between the two non-adjacent wells. 

4 Conclusion 

While tunnelling in one dimensional time-independent systems can hardly of- 
fers some surprises, it remains sufficiently rich to provide a warming up for 
dealing with much more elaborate situations that occurs in higher dimensions 
(or when adding an explicit time dependency). In the latter cases, tunnelling 
frequencies or rates generically fiuctuate by several order of magnitudes; under- 
standing these fiuctuations, that can be seen as an intricate overlap of resonan- 
ces — still qualitatively very different from what occurs in one-dimension — , is 
the aim of a vivid field of actual researches both theoretically and experimen- 
tally (Keshavamurthy & Schlagheck, 2011) and many issues still remain unclear. 
The two three-well models we have studied in this paper present the simplest 
bounded situation where tunnelling deviates from the usual purely exponential 
behaviour and is indeed drastically enhanced near resonances. 



El 
E3 



E-V2S 



(62) 
(63) 
(64) 



E 



E + V2S. 
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